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EXECUTIVE  SUMMARY 
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Reference  1.  The  objective  here  was  to  develop  an  alternative  method  based  on 
a  state  space  approach.  This  has  been  done  and  is  described  herein.  Example 
problems  illustrating  the  application  of  the  theory  are  also  included,  in  the 
Appendix,  with  simulation  results  to  verify  the  theory.  The  method  developed 
here  retains  the  drawback  of  earlier  methods  [1]  in  that  it  is  not  always 
possible  to  stabilize  the  non-critical  variables.  This  is  demonstrated  in 
Example  1  of  Appendix.  This  lack  of  stabilization  is,  however,  easier  to  spot 
than  in  the  older  method. 
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I.  INTRODUCTION 


Many  real-world  control  problems  require  an  ability  to  maintain  adequate 
control  in  the  face  of  disturbances.  Often,  this  can  be  handled  by  applying 
the  complete  cancellation  condition  of  Disturbance  Accommodating  Control  (DAC) 
theory.  The  complete  cancellation  condition  allows  us  to  cancel  the  distur¬ 
bance's  effect  on  all  the  state  variables.  This  condition  can  not  always  be 
met;  thus,  a  method  is  needed  to  cancel  the  disturbance's  effect  on  certain 
critical  variables.  A  method  for  doing  this  is  described  in  Reference  l. 

But,  as  it  is  sometimes  difficult  to  apply,  an  alternate  method  was  sought. 
This  new  method  I3  developed  and  described  herein. 

II.  ALTERNATE  METHOD  OF  DISTURBANCE  ABSORPTION  FOR  CRITICAL  VARIABLES 

In  this  report,  we  are  concerned  with  finite-dimensional  controlled 
linear  systems  governed  by  differential  equations  of  the  forms: 


X  *  A*X  +  B*U  +  F*W  (I) 


where  X  »  (Xi,...,Xn)  is  the  state  vector  for  the  system,  U  *  (U 1 . Ur)  is 

the  control  input,  and  W  *  (Wi,...,Wp)  is  the  disturbance  input.  The  matrices 
A,  B,  and  F  are  assumed  to  be  known  and  constant.  The  class  of  disturbances 
to  be  considered  consists  of  waveform  structured  disturbances  [2]  which  can  be 
modeled  by: 


W(t)  - 

H*Z 

(2) 

z 

D*Z 

(3) 

where  the  vector  Z  -  <Zi , . . . ,Zy)  ls  the  state  of  the  disturbance  W.  The 
matrices  H  and  D  are  known  and  constant.  The  disturbance  is  not  directly 
measurable,  in  general,  but  can  be  estimated  using  the  techniques  available  in 
DAC  theory  described  in  the  references. 

The  problem  is  to  stabilize  the  state  vector  X  to  a  linear  subspace  such 
that  the  critical  state  variables  Xc  -  (Xcl , . . . ,Xcm) ,  m<n,  are  asymptotically 
stable  to  the  null  solution  Xc(t)  -  0  with  the  remaining  (non-crltlcal)  state 
variables  remaining  bound  in  the  subspace.  Assume  that  the  complete  cancel¬ 
lation  condition  of  DAC  theory  (B*r  +  F*H  =  0)  is  not  possible  for  any  choice 
of  r.  Thus,  the  general  disturbance  cancellation  theory  cannot  be  applied  to 
the  problem. 


Begin  by  partitioning  the  system  given  in  Equation  1  into  critical  and 
non-crttical  variables  as  shown  below: 


• 

Xc 

All 

Al2 

Xc 

+ 

Bl 

U 

+ 

Fl 

Xnc 

A21 

A22 

Xnc 

B  2 

F  2 

Since  it  is  desired  that  the  critical  state  variables  (Xc)  be  stabilized  to 
zero,  the  control  needs  to  be  designed  to  make  the  critical  state  variables 
asymptotically  stable  while  canceling  any  effects  of  the  disturbance  or  non' 
critical  states  on  the  Xc.  Thus,  U  could  be  designed  as: 


U  «  Up  +  Ud  «  Kp*Xc  +  ri*Z  +  f2*Xnc 


(5) 


where  Up  *  Kp*Xc  is  responsible  for  stabilizing  the  critical  variables  to  zero 
while  Ud  *  n*Z  +  T2*Xnc  cancels  the  effects  of  the  disturbance  and  non- 
crttical  variables  on  the  Xc.  This  requires: 


Bl*(n*Z  +  F2*Xnc)  +  Al2*Xnc  +  Fi*H*Z  -  0  (6) 


Bl*ri 

-  -Fl*H 

(7) 

Bi*r2 

*  -A12 

(3) 

where  Kp  is  picked  to  stabilize: 


Xc 


(All  +  Bl*Kp)Xc 


(9) 


Note  that  the  non-critlcal  state  variables  are  Ignored,  and  this  is  unwise 
unless  they  are  Inherently  stable,  with  respect  to  the  subspace.  Also,  they 
must  remain  bounded  in  the  subspace.  If  the  non-critical  variables  start  to 
grow  without  bound,  most  real  system  will  soon  run  into  trouble.  Unfortunately, 
there  is  not  much  that  can  be  done  to  stabilize  them  using  the  above  procedure. 
It  may  sometimes  be  possible  to  exploit  the  non-uniqueness  of  the  control  action 
in  stabilizing  the  non-critical  variables.  The  general  form  for  the  control 
gain  has  been  developed  in  Reference  3  and  can  be  combined  with  the  modified 
disturbance  cancellation  condition  In  Equation  7  to  yield  a  more  general  form 
for  the  control  parameters.  The  design  equations  for  the  control  parame¬ 
ters  then  become: 


Bl*r  -  -Fl*H 


(10) 
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K  -  -(C*B)f*C*A*M*M#  +  [I-(C*B)f*C*B]<j,*M# 

+  9 [ I  -  M*M^]  (11) 

where 

C  Is  an  m  x  n  matrix  such  chat  Xc  *  C*X 
M  is  an  n  x  n-m  matrix  of  rank  n-m  spanning  the 
subspace  such  that  C*M  =  0 
(!)f  is  the  Moore-Penrose  generalized  inverse  of  (!) 

*  (M'*M)"^*M'‘  where  M'  is  the  transpose  of  M 
<j>  is  an  r  x  n-m  real  parameter  matrix 
9  is  an  r  x  n  real  parameter  matrix  chosen  to  make  the 
critical  variables  asymptotically  stable 
I  is  the  Identity  matrix 

The  control  action  is  given  by  U  *  K*X  +  r*Z  with  K  and  r  found  from  the  above 

relations.  Note  that  using  the  K  given  in  Equation  11  automatically  cancels 

the  Al2*Xnc  term  thus  doing  away  with  the  need  for  Equation  8.  When  the 

generalized  Inverse  of  C*B  multiplied  with  C*B  is  not  equal  to  I,  we  can 

attempt  Co  specify  $  in  the  second  term  in  Equation  11  to  stabilize  the  non- 
crttlcal  variables. 

When  the  above  procedure  is  examined  closely,  it  becomes  apparent  that  it 
fails  for  Bi  =  0.  This  is  due  to  the  dependence  on  U  to  cancel  the  effects  of 
the  disturbance  and  the  non-crltlcal  variables  on  the  critical  variables.  If 
Bl  is  zero  the  control  cannot  directly  affect  the  critical  variables  and  thus 
the  only  alternative  is  to  maneuver  the  non-crltlcal  state  variables  to  cancel 
the  disturbance  effects  on  the  Xc.  Thus,  it  is  desired  that: 

Ai2*Xnc  +  Fl*H*Z  -0  (12) 

To  accomplish  this,  a  servo  tracking  system  is  designed  where  the  servo  com¬ 
mand  drives  Xnc  to  cancel  the  disturbance  effects  on  the  critical  variables. 
Define  the  servo  tracking  error  as: 

e  -  X  -  Xs  (13) 


with  the  servo  command: 
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Which  when  written  out  becomes: 


ec  *  Aii*ec  +  Ai2*enc  +  Ai2*r*Z  +  Fi*H*Z  (17) 

•  • 
enc  “  A21*ec  +  A22*enc  +  B2*U  -  r*Z 

+  A22*r*Z  +  F2*H*Z  (13) 


Now  split  the  coutrol  action  into  U  *  Up  +  Ud  and  use  Ud  to  cancel  the  distur 
bance  terras.  Cancellation  of  the  disturbance  effects  on  e  requires  that: 


Ai2*r 

and 

B2*Ud 

Thus ,  choose  Ud 
B2*T2 

and 

B2*T3 


-Fl*H 


-A22*r*Z  -  F2*H*Z  +  r*Z 
T2*z  +  T3*z  such  that: 
_A22*r  -  F2*H 


(19) 


(20) 


(21) 


(22) 


where  r  must  satisfy  Equation  19.  To  obtain  zero  tracking  error  It  is 
required  that: 


€ 


(A  +  B*K)e 


(23) 


be  asymptotically  stable.  Thus,  design  Up  *  K*e  to  place  the  eigenvalues  of 
A  +  B*K  in  the  left  half  plane. 

Note  that  if  some  of  the  rows  of  32  are  also  zero,  then  they  are  not 
directly  controllable  by  U  and  another  servo  tracker  may  be  needed  to  control 
them  to  Che  desired  values.  This  is  demonstrated  in  Example  3. 

Equations  19,  21,  and  22  can  be  combined  to  yield: 


Al2*B2*T3  — Fl*H 


and 


(24) 


B2*T2  *  ~A22*®2*r3  ”  F2*H 


(23) 
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Replace  r  In  Equation  14  with  r  *  B2*T3.  This  eliminates  the  need  to  find  r* 
Now  apply  the  generalized  Inverse  to  obtain: 

n  *"(A12*B2)f  *Fl*H  (26) 

T2  -*(B2)f  (A22*B2*T3  +  F2*H)  (27) 

Equations  26  and  27,  along  with  satisfactory  stabilization  of  Equation  23, 
are  sufficient  for  a  solution.  They  are  not  necessary  however,  as  can  be 
easily  seen  from  Example  3  where  a  satisfactory  solution  Is  obtained  when  no 
solution  exists  for  26  and  27. 

Note  that  this  procedure  depends  on  the  ability  to  maneuver  the  non- 
crltlcal  variables  as  desired.  This  may  not  be  possible  In  some  cases,  while 
In  others  a  surplus  of  maneuverable  variables  may  be  available  giving  the 
engineer  some  design  freedom. 


III.  DESIGN  PROCEDURE  SUMMARY 


The  design  procedure  Is  summarized  for  the  two  cases  developed  previously. 
First,  partition  the  system  Into  critical  and  non-crltlcal  state  variables  as 
shown : 


Xc 

All 

Al2 

Xc 

+ 

B  1  ‘ 

U  + 

‘  F  1  ‘ 

.  Xnc 

„  A21 

A22  . 

B2  _ 

_  F2  _ 

If  Bi  Is  not  equal  to  zero,  then  apply  the  first  procedure  and  attempt  to 
design  the  control  action  as  U  *  K*X  +  r*Z  where  K  and  r  must  satisfy: 

Bl*r  -  -Fl*H 

and 

K  -  -(C*B)f  *C*A*M*M#  +  [I-(C*B)f  *C*BJ<f,*M# 

+  0 [ I  -  M*M#] 

where  i 


C  Is  an  m  x  n  matrix  such  that  Xc  »  C*X 
M  Is  an  n  x  n-ra  matrix  of  rank  n-ra  spanning  the 
subspace  such  that  C*M  =  0 

(!)f  Is  the  Moore-Penrose  generalized  Inverse  of  (!) 

M#  »  (M‘,*M)“^  *M"  where  M'  Is  the  transpose  of  M 
$  is  an  r  x  n-m  real  parameter  matrix 
9  Is  an  r  x  n  real  parameter  matrix  chosen  to  make  the 
critical  variables  asymptotically  stable 
I  Is  the  identity  matrix 
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Note,  this  procedure  Is  not  always  successful  even  though  K  and  r  are  obtained 
meeting  the  above  criteria.  This  method  depends  on  the  non-critical  state 
variables  remaining  stable  and  bounded  in  the  subspace.  This  may  occur 
naturally  or  when  C*B  pre-mult Iplled  with  its  generalized  inverse  Is  not  equal 
to  I,  we  may  be  able  to  specify  41  above  such  that  It  aids  in  stabilization  of 
the  non-crttical  variables. 

For  Bl  =  0,  use  the  second  procedure.  Design  f2,  T3  to  satisfy: 
n  »-(Al2*B2)f*Fl*H 
T2  «  -(B2 )f(A22*B2*T3  +  F2*H) 
and  Up  *  K*e  to  stabilize 
• 

e  -  (A  +  B*K)e 

The  final  control  is  then  given  by: 


K*(X  -  Xs)  +  T2*Z  +  T3*Z 


where 


Xs 


0 

32*T3 


Note  that  meeting  these  conditions  are  sufficient  for  a  solution,  but  they  ars 
not  necessary.  In  particular.  If  the  B2  matrix  has  some  zero  rows,  the  servo 
tracking  Idea  can  be  repeated  to  maneuver  these  indirectly  controlled  states 
and  possibly  meet  our  design  criteria. 


IV.  CONCLUSIONS/REMARKS 


A  new  design  procedure  for  critical  variable  disturbance  absorption 
controllers  has  been  presented  here.  This  procedure  consists  of  two  basic 
methods  as  described  above  and  is  demonstrated  in  the  included  examples. 
Simulation  of  the  examples  was  done  to  verify  the  design  and  no  problems  were 
noted.  This  method  retains  the  drawback  of  the  earlier  method  (Reference  1) 
in  that  it  is  not  always  possible  to  stabilize  the  non-critical  variables. 
This  is  demonstrated  in  Example  1.  This  lack  of  stabilization  is,  however, 
easier  to  spot  than  in  the  older  method.  The  second  procedure  depends  on  the 
ability  to  maneuver  the  non-critical  variables  which  is  very  application 
dependent  and  may  or  may  not  be  possible. 

It  is  hoped  that  further  researchers  will  carry  on  the  development  and 
refine  the  method  still  more  with  hopefully  the  addition  of  necessary  and 
sufficient  conditions  for  a  successful  application  of  the  design. 
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APPENDIX 


EXAMPLES 


EXAMPLE  #1 


Objective:  Design  the  control  input  (U)  to  stabilize  the  given  system 

where  Xx  is  the  critical  variable  X,  and  Xa  is  the  noncritical  variable  Xnc> 


1  1 


0  a 


U  +  W  5  X.  =  ci  0]X 

1  1 


( System ) 


W  is  a  piecewise  constant  disturbance  and  a  is  a  constant.  Assume: 
W  =  Z  ==>  H=l 


Z  =  <7(  t)  ; 


with  o(t)  zero  almost  everywhere. 


Partition  the  system: 


i  :  i  x 


o  :  a  x„ 


B1  is  not  zero,  so  try  the  control  action  U  *  K*X  +  r*Z  with  r  and  K  given  by 
equations  10  and  11  respectively. 

Bi*r  =  — fv*h  **>  i*r  *  -0*1  «>  r  =  o. 


where 


K  *  -(C*B) '*C*A*M*M"  +  C I-(C*B) '*C*B]«*M*  ♦  ©Cl  - 


M  =  o  5  C  *  Cl  0] 
1 

. 


thus  0 

li«  =  =  (CO  1]  1  )-»C0  1]  =  CO  1] 


< CB ) *  =  (Cl  0)  1  )'  a  1 


K  »  -Cl  1]  0  0 

0  1 


+  Cl  -  1 ]«*M* 


*  C9i  9=]  1  0 

0  0 


k  =«  co  -l]  +  ce*  o]  =  ce,.  -i] 

The  control  action  is  then: 

U  =*  0x*Xt  -  Xs 

Substituting  U  into  the  original  system  yields: 


xt  *  xt  *  Xa  *  8JX1  -  Xa  =  Xt«(l  +  9,,) 


X*  =  a*Xs  *  e^X,.  -  Xa  +  2 

Now  design  ©,.  to  maka  Xx  asymptotic*  1  1  y  stable.  Pick  ©*  =  -2,  then 
U  =  -2* X x  -  X2 
and  the  system  becomes: 
x  1  *  -  X  i 

X*  *  (a  -  l)Xa  -  2*Xx  +  Z 

Note:  The  stability  of  Xs  depends  on  the  value  of  a.  Simulation 

results  are  shown  for: 

a  *  -1  »>  X3  is  stable 
a  =  2  =** >  X3  is  unstable 


A- 2 


0.00 


EXAMPLE  #2 


Objective:  Design  the  control  input  (U)  to  stabilize  the  given  system 

where  X*  is  the  critical  variable  X«.  Xa  and  X3  are  the  noncritical  variables 
making  up  the  vector  X„«. 


0  1 

0  0 

-1  1 


1 

0  W  5  X.  *  [1  0  0]X 

0 

(System) 


w  is  a  piecewise  constant  disturbance.  Assume: 
W  *  Z  «>  H  =  1 


and 

Z  *  o ( t ) ;  with  a(t)  zero  almost  everywhere. 

Salulian  l 

Partition  the  system: 


X. 

o  :  l 

0 

X« 

0 

1 

# 

3 

o  :  o 

l 

4* 

1 

u  + 

0 

-l  :  l 

0 

x„« 

1 

0 

- 

• 

. 

■  • 

. 

. 

Bi  3  0  for  this  system,  thus  we  use  the  second  procedure  and  design 


ud  •  ra*z  +  ra*z 


where  rs  end  r3  are  given  by  equations  26  and  27. 
T*  *  -(Ata*Ba) '«Pi*H  »-Fa*H  *  -1 

ra  »  -(Ba)'(Aaa*Ba*ra  «■  fs*h)  -  -i/2ci  i] 

Note  that: 

<Ba)'  *  (Ba- Ba)-»-Ba'  *  1/2C 1  13 

The  control  action  U«  is: 

U.  *  Z  -  Z 


Now  design  U»  »  K«  to  stabilize  equation  23: 

«  *  <A  ♦  B*K ) « 

Use  the  third  order  ITAE  response  as  the  model  equation. 

£■»  ♦  1 . 73«£3  ♦  2 . 13*£  ♦  1  -  0 
This  yields  X  *  [-1.4  -0.33  -1.4]  and  since 

U.  «  K*(X  -  X.) 


A- 5 


where 


The  total  control  action  is  then: 

U  *  U„  *  U.  *  -1 . 4X  t  -  0.33<Xa  +  Z)  -  1.4(X3  +  Z )  ♦  Z  -  Z 
For  M  piecewise  constant  Z  a  0  almost  everywhere.  Therefore: 

U  =  -1.4X,.  -  0.33Xa  -  I.4X3  -  0.73Z 

Now  design  a  full  order  estimator  to  estimate  Z.  The  composite  system  is  giver 
by : 


X 


0  10  1 
0  0  10 
-110  0 
0  0  0  0 


X  ♦ 


and  the  observer  with  poles  placed  at  -4  is: 


• 

-16  101 

0 

16 

-333  010 

y* 

1 

333 

X  * 

-272  100 
236  0  0  0 

x  * 

1 

0 

u  ♦ 

271 

-236 
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Figure  A-4.  Noncritical 


EXAMPLE  #3 


Objective!  Design  the  control  input  (U)  to  stabilize  the  given  system 
where  XA  is  the  critical  variable  X«.  Xa  and  X3  are  the  noncritical  variables 
making  up  the  vector  X„«. 


. 

0 

1 

0 

0 

1 

X  =* 

0 

0 

1 

X  + 

0 

u  + 

0 

0 

1 

0 

1 

0 

. 

. 

„ 

• 

W  is  a  piecewise  constant  disturbance. 


Assume : 


W  =  Z  =**>  H  *  1 


and 

Z  *  a(t);  with  a( t)  zero  almost  everywhere. 

Solution; 

Partition  the  systems 


— 

w 

r 

■ 

X« 

o  :  i  o 

i 

•  i 

x  1 

0 

1 

• 

s 

o  ;  o  l 

♦ 

0 

u 

0 

0  110 

1 

0 

- 

m 

. 

. 

Bi  s  o  for  this  system,  thus  we  try  the  second  procedure  and  design 

ud  =  ra*z  +  r3*z 


where  r=  and  r3  are  given  by  equations  26  and  27. 
r3  =  -(A13*Ba) 'IFilH 


ra  *  -(Ba) ' ( Aa=*Ba*r3  ♦  Fa*H) 
Note  that; 


A^a*8a  ™  0 


Inspection  of  these  two  equations  reveals  no  solution  is  possible.  Since  B= 
has  a  zero  row,  let's  cry  to  servo  control  Xa  also.  Use  equation  19  to 
determine  r  and  substitute  r  into  the  equation  for  the  servo  tracking  error. 
Design  T  to  satisfy  Ata*r  =  -FVIH 


Cl 

Equation  16 


0]T  *  -1  «>  r 


where  t  is  arbitrary 


is  used  for  the  servo  tracking  error.  Substituting  in  r  yields: 


«, 

+  i 

0 

U  > 

A1=*r  *  p j. *m] 

z  - 

0 

Bs 

ft=_*r  -  f-*hj 

r 

A- 10 


The  control  cannot  completely  cancel  tne  effect  on  the  states,  let's  try 

using  «•*  to  cancel  the  disturbance  acting  on  in  a  manner  similar  to  what 

was  done  with  the  critical  state.  Define: 

0 

|j  =  e  -  s„  with  =  0 

l?*Z  +  a*Z 

Note  that  the  set  point  state  must  now  contain  the  derivative  of  the 
disturbance  since  it  appears  in  the  6  state  equation.  Taking  the  derivative 
yields: 

A  *  «  - 

which  in  matrix  form  is: 


or  when  written  out  becomes: 


Ax  = 

As 

• 

• 

• 

Ms  * 

A -i 

♦  0*Z 

>  a*Z  *  t*Z  - 

z 

. 

•  * 

•  • 

As  = 

As 

+  u  - 

Z  +  T*z  -  |3*Z 

«■  atl 

T,  (3,  and  a  must  be  chosen  to  remove  the  disturbance  on  62  which  is  equal  to 
p2.  Thus  letting  (3  =  t  -  0  removes  the  disturbance  and  setting  a  =  1  removes 
the  derivative  of  the  disturbance.  U«  can  be  used  to  remove  the  remaining 
disturbance  terms  (those  in  the  equation  for  p3).  Noting  that  (3  =  t  =  0  and 
a  *  1,  makes  it  matter  of  inspection  to  determine  Ua  as: 

•• 

u.  =  z  -  z 

Now  we  have  removed  the  disturbance  from  the  system  and  desire  to  stabilize  it 
using  U».  Let: 

U.  *  Kp 

and  determine  K  to  stabilize: 
p  *  (A  ♦  8*K)p 

Use  the  third  order  ITAE  response  as  the  model  equation. 


A- 11 


£  -  +  1.75*£=  +  2 . 15*£  +  1 


0 


Then 

U.  =  -pi  -  3.13*pa  -  1 . 75*  ps 
But  p  =  6  -  e.„t  therefore: 


Mi  =  «i 
Mi  =  *2 

Mi5  =  «3-(3*Z-a*Z  *  -  Z 


So 


-«i  -  3.13**a  -  1.75*«3  ♦  1 . 73* Z 


But  «  =*  X  -  Xmm,  therfore: 

«i  =*  Xi 

«*  »  X3  +■  Z 
*3  =  X3 


which  yields; 


U.  *  -X 1  -  3 . 13*  Xa  -  3 . 1 3* Z  -  1.75*X,  *  1.75*Z 


The  total  control  action  is  then: 

•  •  ♦ 

u  *  U,  ♦  U„  *  -Xi  -  3 . 13Xa  -  1.73X*  -  2 . 13* Z  ♦  1 . 75* Z  -  Z 

•  «e 

For  W  piecewise  constant  Z  s  0  and  Z  s  0  almost  everywhere.  Therefore: 

U  »  -Xi  -  3 . lSXa  -  1.75X3  -  2.15Z 

Now  design  a  full  order  estimator  to  estimate  Z.  The  composite  system  is  giver 
by : 


« 

0  10  1 

0 

Xi 

0  0  10 

0 

X* 

X  = 

0  10  0 

X  + 

1 

U  5  X  = 

Xi 

0  0  0  0 

0 

z 

. 

_ 

. 

and  the  observer  with  poles  placed  at  -4  is: 


-16  101 

0 

16 

-333  010 
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Figure  A-7.  Noncritical  variables. 
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Figure  A-8.  Disturbance . 
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